Nonequilibrium interacting electrons in a ferromagnet 
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Dynamics of the magnetization in ferromagnets is examined in the presence of transport electrons 
allowing the latter to interact. It is found that the existence of inhomogeneities such as domain wall 
(DW) structures, leads to changes that affect the dynamical structure of the equations of motion 
for the magnetization. Only in the limit of uniform magnetizations or sufficiently wide DW's, the 
equations of motion maintain the form they have in the noninteracting case. In this limit, results 
like the spin torques, the Gilbert parameter, and the DW velocities become renormalized. However 
the length scale that defines such a limit depends on the strength of the interaction. It is shown that 
if large ferromagnetic fluctuations exist in the metallic band then the range for which conformity 
with the noninteracting case holds extends to the limit of arbitrarily narrow DW's. 
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Increasing possibilities for spin selective transport of 
electrons in different materials and their promising phys- 
ical and technological consequences have yielded a broad 
variety of research on electronic motion in the presence of 
a ferromagnetic background. Interest ranges from semi- 
conductors to transition metal ferromagnetsii^ii^ii^i^i^ii. 

A question under current attention is how electrons 
flowing in a ferromagnetic media should modify the dy- 
namics of its order parameter, M(r, i). In the recent 
years, a typical route to handle this problem has con- 
sisted of treating conduction and valence electrons dis- 
tinctively, the former assumed to belong to a s-band, 
while the latter being in a filled d-band responsible for 
the local moments that give rise to M. 

One often uses some equation of motion (as, e.g., the 
Landau-Lifschitz-Gilbert (LLG)^) to treat M as a clas- 
sical field, while a kinetic transport equation is used for 
the dynamics of the conduction s-electrons, usually in 
the limit of long length scales where semiclassical equa- 
tions are valid for the electronic spin^. Coupling of these 
equations is achieved by including a s-d channel through 
a Kasuya interaction operator'^?!!^, whose form is famil- 
iar, e.g., in the Kondo physics of dilute magnetic alloys^. 
In transition metal ferromagnets this approach has been 
used in attempts to draw an adequate picture to under- 
stand experimentsi^iiiii^ii^. Recently discussed issues in- 
clude the role of the so called nonadiabatic spin torque 
that the electronic current produces on a domain wall 
(DW), affecting its motion and stability, and the con- 
tribution of the conduction electrons for the damping of 

Underlying the theories used to treat these phenomena 
we may identify two assumptions (among many others). 
First, that the ferromagnetic state is a broken symmetry 
of the interacting electrons within the d-band. Second, 
that the s-band interelectronic interaction (s-s) leads 
only to trivial renormalizations or else that the s-s inter- 
action is sufficiently weak not to produce any measurable 
effects on M. Thus while the interaction is considered to 
be sufficiently important in the d-band to drive its elec- 
trons into a ferromagnetic state, the s-band is taken as 



an essentially noninteracting-electrons' metallic band. 

In this article we relax the second of these assumptions 
and investigate some preliminary general consequences 
that are expected to be seen on the dynamics of the mag- 
netization. As we will see, the interactions may affect 
the dynamical structure of the equations of motion. In 
the limit of uniform magnetizations or sufficiently wide 
DW's, treating the s electrons as noninteracting is justi- 
fied with the proper renormalizations. We show that in 
the particular limit for which large ferromagnetic fluctu- 
ations exist in the s-band, decoupling of the equations 
of motion in a local frame is possible for narrow DW's. 
Consequently, not only the length for which the "wide 
DW" approximations become valid is lowered, but the 
equations of motion acquire the form they have for non- 
interacting s-electrons. 

For these purposes, we take the usual algorithm de- 
scribed in the third paragraph, and modify it to include 
the s-s interactions. This is done by using the Fermi 
liquid theor}«^ i^^i^° and the results should hold provided 
the conduction electrons are normal (no assumption is 
required to whether the interactions are weak or not). 

The transport equation for the s-electrons' distribution 
can be written in a compact form^, 



dtM + {M,n}-^-[n,U] =/[AA], 



(1) 



where all quantities are 2x2 matrices in spin space, 
7Vp(r,t) is the full distribution function for the s- 
electrons, T-ip{r,t) is the effective hamiltonian and /[A/] 
is a collision integral. Curly Poisson and straight com- 
mutator brackets are also implied in the notation. The 
commutator is kept with respect to the spin matrices as 
a consequence of using a local quantization axis. These 
functions must be understood as averages in the state of 
the system, which satisfy quantum mechanical equations 
of motion that reduce to the above semiclassical equation 
in the long wavelength limit^-''. 

Inclusion of the s-d coupling to the order parameter of 
the ferromagnetic band is achived by using an exchange 
interaction besides the Fermi liquid hamilitonain TCfl, 



n = Hfl + Ksd, where Hsd = SJe^a ■ M(r,t)/Af,,S' 
is the projection of the local spin for the d-electrons, 
Ms is the saturation magnetization, and a are the usual 
Pauli matrices'^-. Note that the structure of the trans- 
port equation in spin space carries a commutator even in 
the absence of the s-d exchange coupling, which has long 
been known^, however the intrinsic s-electron's commu- 
tator vanish if these electrons are noninteracting. 

To write equations of motion from the transport equa- 
tion, we trace Eq.([T]) to obtain equations for the local 
spin density of the s-electrons m — tr[(Tj\f] and the spin 
current tensor J^ = tr[cr{dH/dpi)N]. 

Before writing these equations we separate, as in 
Ref.^^, the local spin density in terms of an adiabatic 
component in the direction of M and a transverse (nona- 
diabatic) component 6m, 



Ms 



(5m, 



(2) 



where Uq = ?^o/(l + ^o )i no is the local equilibrium spin 
density, and F" is the usual notation for the standard 
Landau antisymmetric dimensionless parameters. 

The spin current tensor can similarly be split in two 
contributions, 



■^-- 



" Jci 



(3) 



where j'e is the electric current density, e is the electron's 
charge, ^b is the Bohr magneton, and P' = P/{1 + F^) 
with the current's spin polarization given by P. 

We note that the adiabatic responses to M are renor- 
malized by the interaction in exactly the same way as 
the Pauli susceptibility in a paramagnet. This comes 
from the s-d coupling that effectively adds a local contri- 
bution 



SJe. 

2M. 



-M, 



(4) 



to the molecular field hp = h" -I- hp^^ acting on the 
Fermi liquid quasiparticles^S. Formally, the local field 
S = —S'M./Mg plays the role of a "magnetic field" in the 
equations of motion of a uncharged Fermi liquid. 

The equation of motion for the local spin density be- 



dm 



V,J, = —MxSm-—, (5) 

TexMs Tsf 



where sum is implied over repeated indices. Equation 
^ is the continuity equation for the s-spins modified 
by keeping the part of the collision integral that does 
not conserve spin (here we use a relaxation time approx- 
imation with Tsf for scattering processes other than the 
ones mediated by the s-s interaction). The term with 
Tej, = h/SJex comes from the s-d hamiltonian that af- 
fects the molecular field through hp given by Eq.(I3]). 



For the spin current, we obtain. 



dtjai + clViSm 
hM 2m 



r.Ms 



hN{0) 



FS- 



pa 



XJa 



-,(6) 



'^/ 



with cl = (w|/3)(l + F^){1 + F^/3), N{0) the density 
of states at the s-band unpolarized Fermi surface, and 



T'^f' - {l+F,y3)TFL-'+rsf- 



The fact that t's^ ^ Tsf 



results from the internal spin diffusion in the Fermi liq- 
uid, for which spin-current is not a conserved quantity. 
Such a difference exists even in the abscence of interac- 
tion (as is easily seen by putting F^ = 0) although this 
has not been always considered. 

Now, with the usual assumption that the s-elcctron's 
dynamics is much faster than that of M, we look for the 
solution of Eq. ((6|) after the transient part has practically 
vanished out, what occurs after a few relaxation times. 
At these times, the quasi-steady solution consists of the 
precession of the spin current about M, so that the Eq.® 
becomes 



gm X j„,^ + jVj = -DsViSm, 



(7) 



where g = 2T'^j:{FI}-F^/3)/hN{0) measures the strenght 

of the interaction and Ds ~ ^pT'sfi^ + Fq)/3 is the spin 
diffusion coefficient. 

In Eq.([7|) the differences between interacting and non- 
interacting s-electrons become apparent. By setting the 
interaction to zero {F^ — 0) one recovers Fick's relation 
between spin current and spin density. The additional 
term on the left side originates from the precession of the 
s-electrons spin current about the local spin density, the 
same phenomenon that causes the long known Leggett- 
Rice effect in normal "^He at very low temperature s^^'^- . 

Combining Eqs.(l7]) and ^, we obtain 



D'y^Sm 



1 



(5m 



-^m X M 
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dM fieP' 



Ms dt 



eMs 



Ue ■ V)M, 



(8) 



where D^ = £>«/(! 



and 



A 



Ms 



V,(MxV,(5m)+f 2_oj V4M(M-V,;(5m)]. 

(9) 
Equation ([8]) is the central result of this article. Its form 
suggests that since the interaction between s-electrons is 
in general not weak, the use of an independent-electrons's 
approach in the s-band, as it has been common in the 
literature for s-d-coupling models, is not correct in prin- 
ciple. Regardless of the equation of motion satisfied by 
M, one sees by a careful examination of Eqs.® and Q 
that the dynamics obtained for M in the absence of s-s 
interaction will be in general rather different from that 
in its presence. 

To illustrate how the s-s interaction may affect the 
dynamics of M, we shortly indicate how the result given 



by Eq.® changes the widely used LLG equation®, with- 
out questioning its hniits of vahdity^'' , since this specific 
choice does not change the general conclusions we want 
to draw in what follows. 

Interaction of the spin density Sm with M, brings 
an additional torque to the LLG equation, T — 
{l/TexMs)Sm X M. From this torque 5m can be taken 
as a function M and used in Eq.® to give an equation 
for the dynamics of the order parameter M. A more 
detailed calculation of how the dynamics of M changes 
with the s-s interaction will be shown elsewhere, but it 
is clear, by simple inspection of Eqs.® and ([9]) that the 
result obtained will be two quite different equations for 
interacting and non-interacting s-electrons. 

In the limit of vanishing s-s interaction these ex- 
pressions give the known results for the equations of 
motionii. The additional dynamical structure in the in- 
teracting case comes from the term A given in Eq.(l9]). 
Besides this additional contribution, there are also renor- 
malizations of the coefficients in the other terms. We 
emphasize that our result is exact, apart from the usual 
assumptions of Fermi liquid theory and the general hy- 
potheses that are used to validate the s-d model. 

We see that in the limit of strong interactions A may 
introduce changes that lead to different dynamics. As 
a matter of fact, the strongly interacting limit is some- 
what difficult to evaluate, since while n^ = no/{l + Fq) 
becomes negligible in the limit of large F°'s, the exact 
behavior of the product griQ, that governs A is difficult 
to predict without an explicit knowledge of the Landau 
parameters and also the effective mass that enters N{0). 
One can, nonetheless, keep the values of the Landau pa- 
rameters finite and 7^ to study the overall differences 
that A is expected to produce. 

It is immediately noted, however, that A is only differ- 
ent from zero if the gradients of M are not zero, i.e., in 
the presence of inhomogeneities. For uniform or nearly 
uniform magnetizations, one sees that the sole effect of 
the interaction is to renormalize the equations of motion 
through the coefficients P' and Uq. 

In what follows, we will avoid the more complicated 
task of solving the full equation for a finite A (Eq.®), 
and rather concentrate in presenting a precise meaning 
for the limit of nearly uniform magnetization when s-s 
interactions are not zero. 

First we assume that in this limit, the DW's width 
(or, equivalently, the characteristic average width of the 
inhomogeneity) is very large compared to any diffusion 
length of the problem, so that M varies slowly in space. 
Under this assumption, decoupling of Eq.® is possible 
in a local frame in a similar way that it is accomplished 
in the noninteracting case. Taking the z axis of such a 
frame parallel to M yields the following equation. 



V-^(5m_ 



1 



where (5m_ = Snix — iSray, f- ~ fx ^ ify, 
1 



/(r) = 



1 + igriQ 






, (11) 



and A^ = D'g{l + ign'Q)/{T^f^ + iT^J-). Equation (fTO]) 
admits a standard Green's function solution. 



(5in_ 



47r7 |r-r'| ■^ ^ ^ 



(12) 



Let the typical scale over which /(r) varies be w. It is 
clear from Eq. (jlip that w is given by the average distance 
over which M varies in space. For DW structures, w is 
given by the average DW's width. If w » |A| then we 
can take /_(r) out of the integral in Eq. p^ . to obtain 
-(5m_/A2 = /_(r), and, from Eq.dTO]), V^Jm = and 
A = 0. 

Equation ^ then aquires the same form as the one 
used for noninteracting s-electrons, and the torque and 
other quantities can be determined by simply replacing 
the renormalized parameters in the results earlier ob- 
tained in the abscence of s-s interactioi*iiii^. 

The condition w » |A| that validates this "wide DW" 
approximation can be explicitly written as 



w » 



^F'sf 



r(l + F-) 



3^1 + {gn'oY 



(13) 



r,,^ +T-^. 



Sra^ =/-(r). 



(10) 



where r^^ — 'sf 

We see that the characteristic length for this limit to be 
set depends on the strength of the interaction. Again, the 
simple limit of very large interactions is difficult to eval- 
uate since one does not in general know how the Landau 
parameters and the effective mass are exactly related. 

However, other strongly-correlated limits are ap- 
proached when the s-electrons are driven towards some 
instability. The general Pomeranchuk instability^ is set 
when F^/{2i -I- 1) ^ — 1+- There are other correlation 
driven instabilities whose onsets could be studied, like 
the metal-insulator transition or the SDW instability^. 

We will focus here on the Stoner instability that is 
approached when Fq -^ ^ 1+- The onset of such an in- 
stability is characterized by the presence of large ferro- 
magnetic fluctuations in the system of s-electrons. If one 
thinks of a "weakly interacting" s-band it may appear un- 
physical to assume the s-electrons are in such a regime. 
However, we just showed that no physical "weakly inter- 
acting" scenario exists in a non-trivial solution (-F^'s not 
all equal to zero). The electrons are rigorously undistin- 
guishable, s-d models separate them in two groups with 
hopes that this will be an effective descripition. It is 
quite compelling then, that the strong ffuctuations that 
are present in the d-electrons' part of the wave function, 
will be also shown in the s-electrons' part. Such a regime 
is analogous to the so-called "nearly ferromanetic Fermi 
liquid", a regime that, e.g, normal '^He undergoes as its 
temperature is lowered at room pressure^S,. 



In the presence of large ferromagnetic fluctuations, 
Eq. p^ becomes 



w » 



3gno 



(1 + Fo«)^0, 



(14) 



and the effective length defining the "wide DW" limit 
is arbitrarily lowered. It is clear, however, that a lower 
limit exists for the usual approximations to be valid (for 
example, the semiclassical equation of motion is valid for 
long wavelengths), however these scales are usually much 
less than the typical diffusion paths that relate to |A| in 
the limit of vanishing s-s interaction. 

As examples, the renormalizations of the gyromagnetic 
ratio and Gilbert parameter in the LLG equation in the 
regime of strong ferromagnetic fluctuations, are given by 



7 = ^7, and a ^ ra -\ . 



(15) 



where r = M,(l + F^){1 + Tl/T^f)/no. 

In closing, we considered the effect of the interelec- 
tronic (s-s) interaction in a s-band coupled to a ferro- 
magnetic d band, and found that structural changes in 
the dynamics of the order parameter M may occur. 

The theory correctly reproduces the nointeracting limit 
that has been used as the standard approach for electrons 
flowing in ferromagnetic metals, which include transition 
metals and their alloys. We saw that, in the general 



case, the only safe assumption to boldly ignore the s-s 
interaction would be that its effective amplitude is too 
small, which does not appear to be a plausible demand 
for these systems. 

In the particular regime of large ferromagnetic fluctu- 
ations in the s-band, the inclusion of the s-s interaction 
is shown to produce two main effects: 1 - rescaling the 
equations of motion, and 2 - extending the "wide inho- 
mogeneity" limit to arbitrarily narrow domain walls. 

This latter effect may be physiscally understood by 
noticing that in such a regime of large ferromagnetic 
fluctuations the spin density (m) in the s band becomes 
highly susceptible to changes in the magnetization (M) 
(analogously to what happens in a usual paramagnet at 
the onset of a ferromagnetic instability) . The closer the 
instability is, the stronger the tendency of m to track M, 
so that the length scales within which M changes must 
be narrower in order for the dynamics of m to be affected 
by them. This is appreciated only if the interactions are 
properly accounted for. 

The fact that early approaches that do not consider the 
s-s channel work reasonably well in explaining some ob- 
served properties ^^"iSi, may be suggestive that the s-band 
electrons in the theoretical description of these materials 
must be considered as a strongly correlated state with 
large ferromagnetic fluctuations. 
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